Symmetry of Magnetic Quantum Tunneling in Single Molecule Magnet Mni 2 -acetate 
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The symmetry of magnetic quantum tunneling has been studied in the prototype single molecule 
magnet Mni2-acetate using a micro-Hall effect magnetometer and superconducting high field vector 
magnet system. An average crystal fourfold symmetry is shown to be due to local molecular envi- 
ronments of twofold symmetry that are rotated by 90° with respect to one another, confirming that 
disorder which lowers the molecule symmetry is at important to magnetic quantum tunneling. We 
have studied a subset of these lower (twofold) site symmetry molecules and present evidence for a 
Berry phase effect consistent with a local twofold symmetry. 

PACS numbers: 75.45.+j, 75.50.Tt, 75.60.Lr 
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Magnetic quantum tunneling (MQT) in molecules is 
a very active area of research, starting with the ini- 
tial observations of resonant quantum tunneling in single 
molecule mag net (SMM) Mnia-acetate [H0. Quantum 
tunneling enables the creation of superposition states 
necessary for quantum computing and is also a mech- 
anism for the loss of information in classical magnetic 
information storage. Condensed matter systems like 
these 0, 0, with discrete energy levels coupled to the 
environment in the solid, are also important to un- 
derstanding the border between quantum and classi- 
cal physics, and the decoherence of quantum systems 
The symmetry of the molecule is the first step in 
the understanding its energy level structure and quan- 
tum tunneling. For example, the symmetry determines 
the allowed terms in the spin Hamiltonian and, in par- 
ticular, the form of the transverse interactions that 
lead to tunneling. However, in Mni2 as well as other 
SMMs MQT does not reflect the molecule site symmetry 
and the origin of tunneling has been an o pen question 
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Mni 2 -acetate [Mni 2 Oi2(CH 3 COO)i 6 (H20)4]- 
2CH3COOH-4H2O consists of a core of twelve manganese 
atoms with a ground state spin of 10. The spin Hamil- 
tonian of an ideal (S4-symmetry) Mni2-acetate molecule 
is: 



H 



-DS 2 z -BS 4 z +C(Si + Si)-g[i B Il-S. (1) 



The first two terms represent the uniaxial magnetic 
anisotropy of the molecule (D > and B > 0). The 
parameters, D = 0.548 K, B = 1.1 x 10~ 3 K and 
C = 3xlO _5 K, have been determined by inelastic neu- 
tron spectroscopy ^lj an d EPR [2 experiments. There 
are thus 25+1 = 21 allowed projections of the spin on 
the z-axis, with up and down projections along the z- 
axis separated by an anisotropy barrier of ~ DS 2 + BS 4 
(~ 66 K). An applied longitudinal magnetic field, H z , 
shifts the energy levels favoring the projections of the 
magnetization parallel to the field direction. There are 
values of the z-axis field (resonance fields) for which the 
levels m and m! with antiparallel projections on the z- 



axis are nearly degenerate, ~ kD/gfiB, k — m + m! , 
(&0.44T) and transitions across the energy barrier can 
occur. Interactions that break the axial symmetry and 
mix the levels m and ml lead to an energy difference 
between symmetric and antisymmetric linear combina- 
tions of spin-projections known as the tunnel splitting. 
The lowest order transverse term allowed by the tetrag- 
onal symmetry is fourth order (i.e., the third term in 
eq. 1) and would lead to the tunneling selection rule 
to — ml = 4i, with i an integer. 

This is not observed experimentally. In addition, re- 
cent experiments show that there is a distribution of tun- 
neling splittings in Mni2 crystals associated with 
disorder. Two distinct models of disorder have been pro- 
posed. Chudnovsky and Garanin 0] proposed that ran- 
dom line dislocations in a crystal lead, via magnetoelastic 
interactions, to a lower molecule symmetry and a broad 
distribution of tunneling rates. Subsequent magnetic re- 
laxation experiments indeed showed the existence of a 
broad distribution of tunneling rates and were analyzed 
in terms of this model 0, ^J. In contrast, Cornia et 
al. [l5j suggested, based on detailed x-ray analysis, that 
variations in the position of the two hydrogen-bonded 
acetic acid molecules surrounding the Mn i2 clusters lead 
to a discrete set of isomers with lower symmetry than 
tetragonal. Our recent magnetic relaxation experiments 
in a longitudinal (easy axis) magnetic field are consistent 
with solvent disorder However, direct evidence of 
lower symmetry molecular environments in MQT exper- 
iments was lacking. 

In this letter we present studies of the symmetry of the 
magnetic response in deuterated Mni2 single crystals in 
the pure quantum regime (T = 0.6 K), in which relax- 
ation is by MQT without thermal activation 0, [2(| . An 
average fourfold symmetry is unambiguously shown to be 
due to local molecular environments of twofold symmetry 
that are rotated by 90 degrees with respect to one an- 
other, in accord with the model proposed by Cornia |l5j| . 
Further, we have studied a subset of these lower site sym- 
metry molecules and present evidence for a Berry phase 
effect consistent with a local lower (twofold) symmetry. 
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The magnetization component parallel to the axial di- 
rection (z-axis) of a Mn 12 -acetate single crystal was mea- 
sured using a high sensitivity micro-Hall effect magne- 
tometer in a low temperature Helium 3 system [5l|. A 
single crystal was placed with one of its faces parallel to 
the plane of the sensor while a high field superconducting 
vector field magnet was used to apply magnetic fields at 
arbitrary directions with respect to the crystallographic 
axes of the sample. To study MQT rates we sweep the ap- 
plied z-axis field at a constant rate (a — dH z /dt) through 
a resonance and measure the change in sample magne- 
tization, Mtefore — Mafter 22]. The normalized mag- 
netization change (M be f ore - M a f ter )/{M be f ore - M eq ), 
where M eq is the equilibrium magnetization, is the MQT 
probability, P. For a monodisperse system of molecules, 
this probability is related to the quantum splitting of 
the resonance, through the Landau-Zener formula 
Plz = 1 - exp(— 7rA^/2^ a), where u = gn B (2S - k) 
and v a is the energy sweep rate. When there is a dis- 
tribution of quantum splittings in a crystal this situation 
changes. Now the MQT probability depends on the dis- 
tribution of tunnel splittings of the molecules in the ini- 
tial state m = 10 (the metastable state) prior to crossing 
the resonance. Therefore, it is possible to study different 
parts of the distribution (such as molecules with either 
the largest or smallest tunnel splittings) by appropriate 
preparation of the initial magnetization state, as will be 
described below. 

In our first experiments, two different magnetization 
states were studied. In the first case a) M ini ti a i = M s 
and we study relaxation with the field applied along the 
-z direction (M eq = —M s ). The initial state corresponds 
to all molecules in the state m = 10 and the entire distri- 
bution of quantum splittings (100%) contributes to the 
relaxation as the field is increased along the -z direc- 
tion. Note that for the smallest observable resonances 
the molecules with the largest tunnel splittings within 
this distribution will make the dominant contribution. 
In the second case, b) we start with Mi n iu a i — — 0.4M S , 
having allowed 70 % of the molecules to relax prior to the 
experiment by crossing resonance k = 6 (in the absence 
of an applied transverse field). We are thus examining 
30% of the molecules with smallest tunneling splittings. 
After we prepare these magnetic states we sweep the lon- 
gitudinal magnetic field at a constant sweep rate (a = 0.4 
T/min) from T to about 5 T, while a transverse field of 
0.4 T is applied at an angle, 0, with respect to one of the 
faces of the crystal. We repeat this procedure for angles 
from (f> = 0° to <p = 360°. 

In fig. 1 we show the angular dependence of the 
MQT relaxation probability for both initial configura- 
tions. Resonance k = 6 is shown for case (a) and k 
7 for case (b). In both cases the result is a four- fold 
pattern of maxima in the MQT probability at <f> = 60°, 
150°, 240° and 330°. These maxima correspond to an 
enhancement of the relaxation rates when the transverse 
field is applied in these directions. There is also a one- 
fold contribution due to a small misalignment between 
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FIG. 1: MQT relaxation probability versus the angle, <j> , be- 
tween the applied transverse magnetic field and one of the 
faces of the crystal for resonance k = 6 (with the whole crys- 
tal) and k — 7 (for 30% of the molecules with smallest tunnel 
splittings). The solid line represents the effect of misalign- 
ment between the longitudinal field and the easy axis of the 
crystal. 

the applied field and the z-axis of the crystal which is 
represented by a continuous line in the figure. It is clear 
that both initial magnetization states show the same be- 
havior for the probability versus angle of the transverse 
field. Since these initial configurations represent two dif- 
ferent parts of the tunnel splittings distribution (high 
and low ends of the distribution), this indicates that the 
four-fold symmetry of the MQT probability is a property 
of a significant fraction of the molecules in the crystal. 
This four-fold rotation pattern is qualitatively consistent 
with the fourth order transverse anisotropy term in the 
Hamiltonian of eq. (1). Note that this term leads to two 
hard and two medium magnetic axes in the x-y plane. 
For example, for positive C the x and y (<f> = mr/2, with 
n = 0,1,2,3) axes are hard and x — ±y (<j) = (2n + l)7r/2) 
are medium axes. A magnetic field applied parallel to 
a medium magnetic axis produces a larger tunnel split- 
ting than the same field applied along a hard axis |23j. 
However, the value of the parameter C found by EPR 
spectroscopy, ~ 3 x 10 ~ 5 K, would produce a change in 
tunneling probability of only ~ 15%, much smaller than 
the results observed in the lower curve in fig. 1. 

In order to determine the origin of this four-fold sym- 
metry and whether it is intrinsic to the Mni2 molecule 
we have conducted the following experiment. We have 
selected a fraction of molecules in the crystal by apply- 
ing a selection transverse field (STF), Hstf — 0.6 T, 
aligned with one of the maxima in fig. 1, i.e. parallel to a 
medium magnetic axis, while crossing resonance fc, with 
the longitudinal field swept along +z. Then we study 
this fraction of molecules: we measure the magnetization 
versus longitudinal with the field swept in the opposite 
direction, along — z, in the presence of a transverse field 
of 0.3 T applied at different angles, </>, as was done in Fig. 
1. First, the STF was set at ^select = 60° and in sepa- 
rate experiments both 50% (selection with k = 6) and 
10% (selection with k = 5) of molecules with the largest 
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FIG. 2: MQT relaxation probability versus (j> for resonance k 
= 6 for different selections of the initial state of the sample. 
In the left figure the STF was applied along 4> ae i ect = 60°, for 
different initial states (50% and 10% of the biggest splittings 
of the distribution). The right figure shows the result for a 
preparation of an initial state with 50% of the distribution 
with the orthogonal STF , (f> se i ec t = 150°. The left polar 
plot shows the results for the biggest splittings (10% of the 
distribution) for both selection angles. The polar plot on the 
right represents the calculated behavior of the splitting for 
E = +10 mK (thick solid line), E = -10 mK (thick dashed 
line), E — (thin solid line). 

tunnel splitting were studied for this angle of the STF. 
The same procedure was repeated with the STF in a di- 
rection orthogonal to the first case (<fi S eiect — 150°), with 
50% of molecules with the largest tunnel splittings. The 
behavior of the MQT probability versus the angle, <f>, is 
shown in fig. 2 for resonance k = 6 and for the three ini- 
tial magnetization states described above. The results in 
all the cases show a two-fold rotation pattern of the MQT 
probability, in clear contrast with the four-fold rotation 
pattern observed when there was no transverse field used 
to select the initial state of the sample. When the STF 
was applied at 4> select — 60° (left figure), the results show 
only two maxima at 60° and 240°. On the other hand, 
for the selection at 4> se iect = 150° the probability shows 
maxima at 150° and 330°. In the left polar plot in fig. 
2 we show the results for both angle selections, 4> se iect = 
60° (solid circles) 4> S eiect — 150° (open circles), where the 
effect of the field misalignment has been corrected. 

The observation of a two-fold symmetry in the MQT 
relaxation probability, with a 90° phase difference that 
depends on the direction of the STF, is clear evidence 
that molecules in the crystal have lower than fourfold 
symmetry. These observations are in excellent accord 
with the isomer model of Cornia et al. ^5|. There 
are six isomers in this model, four of which have lower 
symmetry than tetragonal (and comprise about 90% of 
the molecules in the crystal) and thus have a second- 
order transverse magnetic anisotropy, E(S^. — ST). This 
anisotropy gives rise to a two-fold rotation symmetry 
for fields in the x-y plane, with maxima in the MQT 
probability separated by 180°. The maxima occur when 
the applied field is aligned with the medium magnetic 
axis (for E < 0, the x-axis and for E > 0, the y-axis). 



Since the crystal has tetragonal symmetry there must 
be equal populations of isomers with opposite signs of 
E. A change in the sign of E rotates the hard and 
medium axes of the molecule by 90 degrees, explaining 
the average four-fold symmetry observed in the crystal 
(fig. 1). Moreover, the isomers found in x-ray diffraction 
have their hard axes in the following directions, 50-60°, 
140-150°, 230-240°, and 320-330° with respect to crys- 
tal faces [24j], as we observe in Fig. 2. Importantly, this 
observation is not consistent with the dislocation model 
|14| , since in this model the distribution of medium axes 
directions is isotropic, excluding very small regions of the 
crystal near the dislocation cores. The polar plot on the 
right in fig. 2 shows the calculated tunnel splittings for 
resonance k = 6 versus (f>. The solid thin line represents 
the angular dependence of the tunnel splitting for E = 0. 
An estimation of the magnitudes of E needed to explain 
the change in the MQT relaxation observed in figs. 1 
and 2 are: E ~ 0.5 mK (fig. 1 upper curve, molecules 
with the smallest tunnel splittings), E ~ 2.5 mK (fig. 1 
lower curve, the whole distribution), E ~ 2.5 mK (fig. 2, 
50% of the largest tunnel splittings with transverse field 
selection) and E ~ 10 mK (fig. 2, 10% with the same 
transverse field selection). Recent EPR experiments sug- 
gest similar values of E 19] . 

By applying a STF in a given direction we can select 
molecules with a particular direction of their medium 
axes. We have thus selected molecules with their medium 
axis along <fi se iect = 60°, and 10% and 50% of the distri- 
bution of such molecules, for further investigation (as in 
fig. 2 on left). We have studied the MQT relaxation as a 
function of the magnitude of the transverse field applied 
along both the hard (cf> = 150°) and medium (<f> = 60°) 
axes of such molecules. The results for resonances k = 5, 
6, and 7 are shown in fig. 3. The probability is plot- 
ted on a logarithmic scale versus the transverse mag- 
netic field. For a given resonance, the difference between 
these orientations of the transverse field is clearly larger 
in the experiment with 10% of the molecules (biggest 
splittings) than with 50%. This indicates that the aver- 
age value of the second order transverse anisotropy in- 
creases as we select the largest tunnel splittings in the 
distribution. An estimation gives average E values of 
£(50%) ~ 2.5ra£, £(10%) - lQmK. These are the 
same values as those extracted from the change of the 
probability in the rotation experiments. One significant 
difference in the transverse field dependence of the prob- 
ability can be observed through the comparison of the 
graphics in fig. 3. In the results for 50% of the distri- 
bution, the probability increases exponentially for both 
transverse field orientations. However, in the experiment 
with 10% of the molecules the MQT relaxation prob- 
ability deviates significantly from exponential behavior, 
showing the largest deviations at fields H p (k = 5) ~ 0.45 
T, H p (k = 6) ~ 0.3 T and H p {k = 7) ~ 0.35 T. This is 
reminiscent of the Berry phase observed in the SMM Fe§ 
|22^ . The modeling that we have carried out using the 
values of E extracted from our data give minima in the 
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FIG. 3: Transverse field dependence of the MQT relaxation 
probability for resonances k — 5, 6 and 7. The transverse field 
is applied along the hard axis (solid symbols) and along the 
medium axis (open symbols). 

tunneling probability for H m i n (k = 5) ~ and 0.5 T, 
H mm (k = 6) ~ 0.25 T and H min (k = 7) - and 0.6 T. 
Some of these minima occur where we observe plateaus 
in fig. 3. The distribution of E values in the experiment 
probably smoothes the minima leading to the plateaus we 
observe. On the other hand, the measured MQT proba- 
bility does not show parity effects for the resonances at 
zero transverse field. In the absence of transverse fields, 
the symmetry breaking terms in the Hamiltonian are due 
to the transverse anisotropy of the molecules. Fourth 



order transverse anisotropy only allows tunneling transi- 
tions for resonances k that are a multiple of 4 (k — At), 
while second order anisotropy only allows transitions for 
k a multiple of 2 (k — 2i). The observation of tunneling 
relaxation at odd resonances must be due to the presence 
of local transverse fields. In the Cornia model this is ex- 
plained naturally as small tilts of the molecular axis of 
the isomers. Other possibilities include transverse dipo- 
lar and nuclear fields or dislocations. 

In summary, the average four-fold symmetry of MQT 
in Mni2-acetate is due to equal populations of isomers 
with biaxial anisotropy that have orthogonal medium 
magnetic axes. Measurements of the tunneling prob- 
ability versus transverse field show evidence for Berry 
phase effects probably due to a combination of second 
and fourth order transverse magnetic anisotropies. From 
a broader perspective, these results illustrate how subtle 
changes in molecule environment can modulate magnetic 
anisotropy and magnetic quantum tunneling. This sensi- 
tivity to local chemical environments is likely to be gen- 
eral to magnetism in molecules and ultimately useful in 
controlling their quantum properties. 
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